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Abstract: We study the spectrum of the gravitational modes in Minkowski space- 
time due to a 6-dimensional warped deformed conifold, i.e., a warped throat, in 
superstring theory. After identifying the zero mode as the usual 4D graviton, we 
present the KK spectrum as well as other excitation modes. Gluing the throat to the 
bulk (a realistic scenario), we see that the graviton has a rather uniform probability 
distribution everywhere while a KK mode is peaked in the throat, as expected. Due 
to the suppressed measure of the throat in the wave function normalization, we hnd 
that a KK mode’s probability in the bulk can be comparable to that of the graviton 
mode. We also present the tunneling probabilities of a KK mode from the inflation¬ 
ary throat to the bulk and to another throat. Due to resonance effect, the latter may 
not be suppressed as natively expected. Implication of this property to reheating 
after brane inflation is discussed. 
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1 . Introduction 

In flux compactihcation in 10-dimensional superstring theory |]I|, such as a KKLT 
vacuum ||^, all moduli of the 6 dimensional compact manifold are stabilized. Pre¬ 
sumably, our today’s universe is described by such a KKLT-like vacuum. Such a 
manifold typically includes warped throats inside the bulk of the manifold. For a 
4-dimensional Minkowski spacetime, a massless 4D graviton is always present as a 
zero mode. In addition, there are excitation modes including Kaluza-Klein (KK) 
modes. It is possible that we (i.e., the standard model particles) live inside a stack 
of D3-branes sitting at the bottom of such a throat. So it will be interesting to 
know the spectrum of the closed string modes in the throat. The spectrum and their 
respective wave functions will reveal the structure of the throat as well as the bulk. 

Another reason in understanding this spectrum and their respective wave func¬ 
tions has to do with the reheating after inflation. It is generally believed that inflation 
explains the origin of our universe. A concrete realization of the inflationary scenario 
in superstring theory is the D-D-brane inflationary scenario [^, ^ . In the KKLMMT 
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scenario 0, where the 6 compactihed dimensions are dynamically stabilized, brane 
inflation is realized by introdncing an extra D3-D3-brane pair in a KKLT-like vac- 
unm. The D3-brane is sitting at the bottom of a Klebanov-Strassler (KS) type throat 
[^], as the mobile D3-brane moves towards it. This describes the inflationary epoch, 
where the vacnnm energy driving inflation comes from the extra D3-D3-brane ten¬ 
sion and a possible conformal conpling term. In this scenario, inflation ends when 
the D3-D3-branes collide and annihilate, initiating the hot big bang epoch. If this 
annihilation happens in a throat (A throat) other than the standard model throat (S 
throat), one likes to ask if the resnlting reheating is compatible with big bang nncle- 
osynthesis. This issne has been studied recently i I, I, 0- A key factor involves 
the tunneling rate, i.e., the probability amplitude of a KK mode (which is produced 
in the A throat) in the bulk and the S throat. Here we examine in more details the 
KK wave functions. 

We start with the KS throat. Given the metric, it is quite straightforward to 
solve the Schrodinger-like equation to obtain the KK spectrum of the lOD graviton. 
We shall assume 4D Minkowski spacetime to be a good approximation to the KKLT- 
like vacuum. We And a zero mode to be identified as the massless 4D graviton, while 
the discrete KK mass spectrum is given by solving a Schrodinger-like equation. 
In the S'^ X approximation around the bottom of the throat, the mass eigenvalues 
may be classified by the corresponding quantum numbers (j — 1) for S'^ and I for 
S^. It is convenient to introduce a dimensionless E^n, A'm ~ EmhA/h where < I 
is the warp factor at the bottom of the A throat and Ig is the string scale. In a 
WKB approximation, which is accurate to around a few percent, we And that 
for m = {n, j — 1,1) (where j — 1 and I are the x S‘^ quantum numbers) is given 
by 



— — Q arctan 



1.03(n-i + ^) 


= 1.4 {f - 1) -K 3.4/(/ + 1) + .65 . 


( 1 . 1 ) 


and the numerical values are chosen to get a good fit to the actual spectrum. The 
lowest is then Ei = 2.02, which, for slow-roll KKLMMT brane inflation |^, takes 
value of order E^ ~ 10“'^Mp. 

Next we like to look at the localization of the KK wavefunctions and their tun¬ 
neling properties. We may restrict ourselves to the S-waves. For j — 1 = Z = 0, we 
have 


-Sm=(n,0,0) — En — 1.04 U -|- 1.01 (1-2) 

We note that both the eigenvalues and eigenvectors are well approximated by Bessel 
functions. The eigenfunctions are given by Bessel functions Ji, of order u ~ 2.44, but 
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with an unusual argument, while the eigenvalues are related to, not surprisingly, the 
zeros of Jy. 

Next we consider the more realistic scenario where the manifold is made of a 
bulk with a number of throats. Treating the bulk crudely, continuity requires the 
matching of each wavefunction and its derivative at the edge of the throat. Now, 
an excitation eigenfunction becomes a linear combination of Jy and Yy in the throat 
and is approximated by sinusoidal functions in the bulk. The wave function of the 
nth KK mode is peaked at the bottom of the throat, as expected. The ratio of the 
wave function T„(rc) at the edge of the throat (where the warp factor /i^(rc) — 1, 
so T„(rc) is comparable to the wave function in the bulk) to that at the bottom of 
the throat 'l/n(0) with warp factor Ha = ^^(0) is given by 

' El {{n + l)hAr (1-3) 


d'n(rc) 

^n(O) 


where the KK mass En — EnhA ~ (n + 1 )/ia (in string units). Another interesting 
quantity is the ratio of \k„(bulk) ~ 'h„(rc) to that of the graviton mode Tg(bulk) at 
the edge of the throat. Since the graviton is everywhere, this provides a measure of 
how much the nth KK mode can get out of the throat when the time scale is longer 
than the tunneling time but shorter than the KK mode lifetime. We hnd that there 
are two limiting situations. In the case where the manifold is made up of a set of 
throats only, 

T„(bulk) 

'hg(bulk) 

This suggests that the tunneling from one throat to another is very much suppressed. 
However, in the case where the manifold has a reasonable bulk (that is, a region where 
there is little or no warping), then 


Et 


((n+ l)hAf 


(1.4) 


T„(bulk) 

'hg(bulk) 


(1.5) 


That is, the average KK wavefunction in the bulk is essentially the same as the 
average graviton wavefunction. This result is due to the fact that the warped metric 
weights the throat signihcantly less than the bulk in the measurement of the volume 
of the compactihed manifold, so the wave function normalization pushes the wave 
function to the bulk. 

On the other hand, the relevant tunneling probability for reheating may be the 
transmission coefficient, which measures the tunneling from the A throat to the bulk 
(and to another throat via the bulk), where there is only outgoing wave (i.e., away 
from the barrier) in the hnal stage. It is important to decide which of the estimates is 
relevant for reheating. Because the compactihed manifold has hnite size, the above 
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ratio estimates may be more relevant if the KK mode decays slowly so its final 
reflection component shonld be inclnded. If it decays fast enongh before reflection, 
then the transmission coefficient shonld be more relevant. In this case, we hnd 

T( A ^ bnlk) ^ ^ ~ ~ ((n + 1) HaT (1.6) 

where 0^/2 is the standard WKB integral of the (imaginary) momentnm across the 
barrier. The transmission coefficient from the A throat to the S throat via the bnlk 
is given by 

T{A ^S)=4 (^(©i + ^)2 cos^ L + 4 sin^ (1.7) 

where En{A) = (n + l)h^ matches En'{S) ~ (n' + 1)^5 in the S throat. In a realistic 
model, the KK modes have widths and this matching is not hard if the warp factor 
of the S throat hs h a- In the case where the manifold is made np of a set of 
throats only (mimicking the two throats model in the RS scenario as in [|^, ^, Q), so 
L ~ 0, we hnd that T is extremely small, 

T(A ^ A) ~ (0^)-" ~ ((n + l)hAY^ (1.8) 


One shonld compare this to that in the Randall-Snndrnm (RS) scenario [0, |T3 


where this ratio is fonnd to go like h\ [0, H, |^. This latter resnlt follows from Ref. []I4| 
where the h\ comes simply from the warped metric. There, the wave fnnction obeys 
the Israel jnnction condition in the infrared region, which is very specihc to the 5- 
dimensional spacetime in the RS model. Here we are dealing with a 6-dimensional 
warped deformed conifold in 10 dimensional spacetime, where the wave fnnction at 
the bottom of the throat mnst be hnite. For Ha ~ 1/5 to 10“^ or smaller, this 
difference is hnge. 

However, for some KK modes, L = {jil + l/2)7r, so cosL = 0. For these KK 
modes, the transmission coefficient (|1.7|) approaches nnity. 


T{A ^ A) ~ 1 


(1.9) 


This is the well-known resonance effect. Based on a simple estimate, we expect that 
some of the KK modes will satisfy this condition. So onr two ways of estimating 
tnnneling of KK modes from one throat to another agree qnalitatively. The end resnlt 
is that the tnnneling depends crncially on the presence/size of the bnlk, bracketing 
the tnnneling amplitnde in the simpler RS scenario. 

The implication to reheating after inflation is interesting. In a realistic Type HB 
compactihed manifold with mnlti-throats, we expect the presence of a bnlk. In this 
scenario, the KK modes are prodnced in the inflationary throat (the A throat) where 
the mobile D3-brane falls into this throat, collides and annihilates with the D3 brane 
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sitting at the bottom of the throat. These KK modes will easily tunnel to another 
throat, say the S throat that contains the standard model branes. These KK modes 
can then decay into the open strings in the S throat. As pointed out in Ref. |]^, 
this will enable the excitation of the standard model modes and reheat the universe 
to start the hot big bang epoch, consistent with big bang nucleosynthesis. The 
differences of our result compared to that in Ref. [0, ^ are the two effects : 

(1) the bound state property in a warped deformed conifold actually suppresses the 
tunneling much more than the 5-dimensional RS scenario; 

(2) the presence of the (un-warped) bulk can provide a huge enhancement of the 
tunneling amplitude. 

These two effects tend to offset each other. In the supergravity approximation, the 
presence of a finite-size bulk is expected, so the second effect may win. In a realistic 
model, under appropriate conditions, we expect the tunneling amplitude can be big 
enough to allow a successful reheating of the universe after brane inflation in a 
multi-throat model. However, the answer may be quite sensitive to the details of the 
model. 


2. Setup and the Graviton 


We are interested in the lOD graviton spectrum. We assume that the background 
is a static geometry with Poincare symmetry in four dimensions. Lacking a deep 
dynamical understanding, a 4D Minkowski metric is achievable only with some hne- 
tuning. There are three kinds of metric fluctuations, labeled by their spin from a 4D 
observer point of view. They are scalar, vector and tensor modes. The scalar modes 
determine the size and shape of the the internal manifold (here Calabi-Yau manifold). 
The scalar modes which control the size of the compactification is known as radions. 
The vector modes are fluctuations of metric in the form g^m, he. metric field with one 
index in the extended four direction and one index in the internal dimension. The 
existence of these modes are model dependent. For example, the Randall-Sundrum 


(RS) model [^| has no massless vector mode due to the Z 2 symmetry. Finally, the 
graviton is the spin 2 fluctuation of the 4D metric, —>• rj^v + h^i,. Beside the zero 

mode, the RS model has a tower of massive KK graviton. These modes are highly 
peaked in the IR region with the mass gap 5m ~ TeV. 

Let us start by considering the following action 


S = M. 


D-2 


d 


D 


X 






topological 


( 2 , 1 ) 


where Mg represents the D-mimensional mass scale. We consider a general form 
of C such that a background solution with a 4D Minkowski metric is possible. C 
contains bulk contributions and it may also contain contributions from local sources. 
For example C in the RS model contains a bulk cosmological constant and localized 
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brane(s) tension(s). In the GKP model, bsed on IIB string theory, C has contribution 
from three form flux G( 3 ), five flux F( 5 ) and axio-dilaton field. It also contains local 
source D3-branes, D7-branes and 0-planes. 

The model admits the following metric 

ds^ = y) dx^ dx"' + gmn{y)dy'^ dy '^. (2.2) 


where the compactihed manifold (with the metric gmn) is stabilize and remains static. 
The Einstein equation is 


Rmn — Tmn — 




(2.3) 


where M and N runs on the entire space-time and T is the trace of the stress energy 
tensor Tmn, given by 


Tmn = 


gMN C — 2 


6C 

6^ 


Suppose C depends only on the internal metric, i.e.. 


6C 


= 0 . 


(2.4) 


(2.5) 


Intuitively this implies that the total contribution from integrating C over the in¬ 
ternal dimensions behaves like a cosmological constant for a 4D observer. Using the 
condition (^), the components of stress energy tensor are 


Tfiy g^vd2 


Tmn gmuR- T Tm 


( 2 . 6 ) 


where T„„ = -liLjiir" . 

With the above form of the stress energy tensor and with the dehnition T = 7)))^, 
the yv components of the Einstein equation obey 




2C + T 
' D-2 


■ gfj.i' 


(2.7) 


It should be emphasized that the above compact form of the Einstein equation is 
true everywhere, in the bulk or at the position of the local sources. At local sources 
like branes, there would be a delta function singularity for C and T such that 
Eq. (|2.7|) is satished. 

With the metric (p.2|) , R^^, is given by 


Ri,u = - ^g^^dmgf.ud'^gpx + ^g'^^dmgupd^gpx . ( 2 . 8 ) 

Here is the 4D Ricci tensor calculated purely from g^^ while Vp_ 4 ) is the 

Laplacian constructed from gmn- 
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For g^v{x,y) = where g^y is Ricci flat {E!'^}{g) = 0), Eg . (12.71) reduces 


to 


(VfD-4,«“ + 5-e^g 

The case we are interested here is 


(2.9) 


9,4 = 




( 2 . 10 ) 


For an extensive analysis of other perturbations, see Ref. K 


In the absence of the perturbation and with the Minkowski background 
A{y) satisfies Eq.( p.9|) . The Minkowski metric implies effective zero 4D cosmological 
constant. In view of the compactiflcation dynamics, this requires a fine tuning. The 
actual fine-tuning required is for a very small positive 4D cosmological constant 
compatible with the observed dark energy. In this case, the Minkowski metric is a 
close enough approximation for our purpose. Next we like to solve for h,. 

Imposing the transverse-traceless gauge conditions 




/i/; = = 0. (2.11) 

we have 

= -^V^4) V - ^Vp_4)(e^^/i^^) - . (2.12) 

where V( 4 ) is constructed purely from g^^. 

Using this expression for and the relation imposed from the background 
solution Eq.(^(^ in Eq (|^) , we obtain 

e"^^Vp_4)(e^^h^^) - e“^^h^^Vp_4)(e^^) -F VJ 4 ) V = V = 0 • (2-13) 


Interestingly enough this is the equation of motion for a D-dimensional massless 
scalar field which was also noticed in [l^. This makes sense because the gravitons 
in D-dimensions are massless spin two particles and since we only excite the /xz/ 
components they behave like scalar fields from the internal dimension point of view. 
This equation is true everywhere in the bulk and the contributions of the local sources 
are included automatically. To obtain the boundary conditions at the positions of 
local sources, one simply integrate passing across the local source. More 

specifically, in a Gauss normal coordinate locally chosen near the source at position 
y = x/o; the boundary condition obtained from Eg. (^.131) is 


d h = 0 


(2.14) 


where is the normal derivative along x* coordinate perpendicular to the source 
where i runs from 1 to the number of transverse dimensions to the source. 
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(2.15) 


The KK spectrum of the graviton is obtained from Eq. (|2.13 ), 

V?B_4)(e''"V) - VV?D-4)(e'^) = 0 

where 


/V = 


(2.16) 


such that Em is the mass measured by a 4D observer. 

It is clear from above equation ( p.l5| ) that there always exists a zero mode with 
solution = constant. This represents the ordinary 4D graviton propagating in 
the four dimension with Poincare symmetry. Since the internal dimensions are com- 
pactihed, the zero mode is normalizable. This in turn ensures that the effective 
four-dimensional gravity strength is hnite. More explicitly, consider the decomposi¬ 
tion of the graviton excitations in the following way 




(2.17) 


The gravitational part is of the action ( p.l|) 
MP-^ 


S„ = — 

^ 2 

Mf-2 

rsj - 

2 

MP-^ 


y X ^/^ R 

j cPx d^~^y \/\gab\ 

[ d^x Y, [ ^''■'2/^^e2^T(m)(|/)4/(mo(|/) (2.18) 


m,m' 


The natural normalization is 

d^-"|/v1^e2^T(m)(2/)^(mo(2/) = M|(5mm' (2.19) 

For the zero mode, Tg, we obtain the usual relation between the gravity strength in 
four dimensions and the fundamental mass scale of the higher dimensional theory. 

After this general discussions we would like to find the KK spectrum of the 
gravitons and other closed string modes in the KS background. To do so, let us hrst 
review the geometry of a warped deformed conifold. 


3. A Throat in the Calabi-Yau Manifold 

A KKLT vacuum involves a Calabi-Yau (CY) manifold with fluxes [^j. Consider F- 
theory compactihed on an elliptic CY 4-fold X. The F-theory 4-fold is a useful way 
of encoding the data of a solution of type IIB string theory; the base manifold M of 
the hbration encodes the IIB geometry, while the variation of the complex structure 















Figure 1: A schematic picture of the Calabi-Yau manifold is presented here. The large 
circle given by dashed line represent the 3-cycle where NS-NS three form is turned 
on. The smaller circle in the throat stands for the 3-cycle where the R-R three form F 3 
is turned on. Also shown are D7-branes wraping 4-cycles. There may exists a number of 
throats like the one shown here. There is a mirror image of the entire picture due to the 
IIB/Z 2 orientifold operation. 


of the elliptic fiber describes the prohle of the IIB axion-dilaton. In such a model, 
one has a tadpole condition on the Euler number x of X 


X(^) 

24 


Xd3 -|- 


2 ^ 10^3 


H?. A iA 


'M 


(3.1) 


Here T 3 is the tension of a D3 brane, is the net number of (D3 - D3) branes 
hlling the noncompact dimensions, and H 3 , F 3 are the 3-form fluxes in the IIB theory 
which arise in the NS and RR sector, respectively. In the absence of flux, it is always 
possible to deform such an F-theory model to a locus in moduli space where it can 
be thought of as an orientifold of a IIB CY compactihcation. So one may use the 
language of IIB orientifolds, with M being the CY 3-fold which is orientifolded. In 
this language, the LHS of Eq. (|3.1|) counts the negative D3-brane charge coming from 
the 03 planes and the induced D3 charge on D7 branes, while the terms on the RHS 
count the net D3 charge from transverse branes and fluxes in the CY manifold. The 
Kahler moduli are stabilized by non-perturbative dynamics. 

As shown in Figure |I], there is a Klebanov-Strassler throat in the manifold [^]. 
The throat is a warped deformed conifolds, which is a non-compact CY 3-fold. This 
throat is glued to the compact CY manifold. A number of such throats can be 
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present in a typical compact CY manifold, although our discussion will focus on a 
specific throat. Such a CY manifold yields a KKLT-like vacuum with a tiny positive 
cosmological constant to explain the observed dark energy. It is expected to be 
metastable, with a lifetime much longer than the age of the universe. 

During inflation, there is an additional D3 - D3-brane pair (without changing 
or any other term in Eq. (|3.1| )). The D3-brane is attracted to a throat, so it is 
expected to sit at the bottom of such a throat. On the other hand, the D3-brane is 
mobile. As the D3-brane moves towards the D3-brane (due to the attractive force 
between them), inflation takes place (due to the extra tensions of this additional 
D3 - D3-brane pair). At the end of inflation, the D3-brane collides and annihilates 
with the D3-brane. As D3-D3 annihilation takes place, the universe is reheated and 
cosmic superstrings are produced. 

A warped throat typically has a non-trivial geometry. The 6-dimensional metric 
around the throat may be written as 

dsg = -I- (3.2) 

where the non-trivial 5-dimensional metric ds^ may take a variety of geometries. 
The best known example is S^, which has N = 4 supersymmetry. However, we 
are interested in Y = 1 supersymmetry, that is, base manifolds that yield CY 3- 
folds. To be specific, we shall focus on the simplest base manifold that has N = 1 
supersymmetry, namely, [1T7|, [I^ . In fact, is the first known Sasaki-Einstein 
metric for N = 1 supersymmetry, which is preserved by its deformation. This conifold 
and its deformation are well-studied |T^, 



3.1 The Conifold 

A cone is defined by the following equation in 

4 

Y “’? = '> (3-3) 

i=l 


Here Eg. (|3.3|) describes a smooth surface apart from the point Wi = 0. The geometry 
around the conifold is studied in Ref.[^. If Wi solves Eq. (|3.3| ), so is uwi for any 
complex u. The base of the cone is a manifold T given by the intersection of the 
space of solutions of Eq.(|3.3|) with a sphere of radius r in = R^, 



= r 


2 


Treating the 4 Wj as a 4-vector w = x -|- zy on with 2 real vectors x and y, then 
the above 2 equations yield 

x-x = rV2, y.y = rV2, x-y = 0 (3.4) 
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The first of these equations defines an with radius r/y/2, while the other 2 equa¬ 
tions define an fiber over the S^. Since all such bundles over are trivial, T has 
the topology of S‘^ x S^. 

We are interested in Ricci-flat metrics on the cone. Consider the ds^ in the 
metric (|^), 

dsl = hahdx°‘dx^ 

is a metric on T, the 6 -dimensional cone admits a Ricci-flat metric if and only if its 
base manifold admits an Einstein metric, that is, Rab{h) = 4:hab- Now let us consider 
a specihc set of ds^ metric: 

2 

ds‘^p,Q = AQ{di/j + Pcos 6*1 d(j)i -|- Q cos 6*2 ^ 02 )^ + A( dO^ + sin^ 6 *j dcfj ) (3.5) 

i=l 

where the integers P, Q are coprime. These are metrics on manifolds which are 

hber bundles over x with P(l) hbers. Now, it turns out that only and 

are S'^ x To be an Einstein metric, requires Aq = 1/9 and ^1 = ^2 = 1/6. 

Furthermore only satishes the Sasaki condition to yield supersymmetry. (A 
Sasaki-Einstein 5-manifold may be defined as an Einstein manifold whose metric 
cone is Ricci-flat and Kahler.) This leads us to study the Sasaki-Einstein manifold 
more detail. 

Let us start with the conifold metric with the base manifold T^’^, 

dsl = d,r‘^ + ^‘^dsipi^i (3.6) 

1 2 ^2 

ds^i.i = -{dtp + cos 6 i dcpi Y + - ( dOf + sin^ 6 ** dcpl ) 

9 ^^ 6 ^^ 

2=1 Z =1 

It can be shown that 

= {SU{2) X SU{2))/U{1) = S^ X SyU{l) 

which has topology of x (with fibered over S^- If and ip 2 are the two 
Euler angles of the two S^s, respectively, then their difference corresponds to P(l) 
while %p = (pi + (p 2 . Since 27i > ipi > 0, the range of ip is An. 


3.2 The Warped Deformed Conifold 


The Klebanov-Strassler throat that we are interested in is actually a warped deformed 
conifold, as illustrated in Figure ^ This warped deformed conifold emerges in the 
presence of fluxes. The R-R flux P 3 wraps the while NS-NS flux P 3 wraps the 
dual 3-cycle B that generates the warped throat, with warp factor h(r). 


1 

47r2/2 



N 


-K, 

KM, 


47r2/2 


F^ = M 


'S3 


hA = e-^^K/3gsM 


(3.7) 
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r=0 



Figure 2: Here is a schematic picture of the conifold (dashed line) and the deformed 
conifold (solid line). The apex is at r = 0. The conifold is deformed at the tip snch 
that r = e is now an S^, where 5^ has shrinked to zero. The dashed circle at constant r 
represents the base of the conifold which is a pgr large r, the base of the deformed 

conifold asymptotically approaches T^d 


Geometrically, the conical singularity of Eq. 
the apex by an 


can be removed by replacing 


'^w^ = e^ ( 3 . 8 ) 

i=l 


where we shall take e to be real and small. The resulting deformed conifold is illus¬ 
trated in Figure ^ and the corresponding metric is non-trivial. It will be convenient 
to work in a diagonal basis of the metric, given by the following basis of 1-forms 

0. 



where 






TT 

-b 

~7r 


2 - e" 

^ ■ V2 

4 

9 = 


V2 




= — sin 01 d(j)i , = dOi , 

= cos Ip sin O 2 d(p 2 — sin xp d 02 , 
e‘^ = sin ip sin 62 dp 2 + cos ip 002 , 

= dip + cos 01 dpi + cos 02 dp 2 


( 3 . 9 ) 


( 3 . 10 ) 
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The metric of the deformed conifold is studied in ^ ^ 

1 


dsl = -e‘'^K(T) 




[SK^t) + cosh^( 2 )[(^ ) + (^ ) ] 

+ sinh2(^)[(^?i)2 + (^?2)^] 


where 


K{t) = 


(sinh(2r) - 2r)^/3 


2^3 sinh 


r 


(3.11) 


(3.12) 


At r —>■ 0, the 5^ (the and components of the metric) shrinks to zero, and we 
are left with 


dnl = i^V3(2/3)V3 (^3)3 ,^4,: 


(3.13) 


which is a spherical S^. Turning on fluxes, the ten-dimensional metric takes the 
following “warped” form 


ds^Q = h dx^dx'^ + ds 


(3.14) 


where dsg is given above in Eq. (|3.11|) . The warp factor h(r) is given by the following 
integral expression 


h(r) = 2^/3 {g.Ma'f /(r), 


(3.15) 


/(r) = 


^^ xcotli x —1 - 2xf^ 

sinh X 


(3.16) 


The asymptotic behavior of K{t) and h(r), which will be important in the following 
calculations, are 


K{t ^ 0) ^ (2/3)1/^ + C>(r2), K{t oo) ^ 

I{t ^ 0) ^ ao + 0(r2), J(r ^ cx)) ^ 3.2-^/3 g-4r/3 

where oq ~ 0.71805. For large r, its relation to r in dsg (|3.2| ) is given by 

3 


^2 _ 4/3 2r/3 


In the limit e —> 0, finite r implies r —> cx), so 


dsl I ^ 


i\2 


9 


(3.17) 


(3.18) 


(3.19) 


2 = 1 
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which is simply the conifold metric (3.6). 

In a realistic model, the throat does not extend indefinitely. At some point 
it should glue smoothly to the bulk of the Calabi-Yau manifold. The bulk of the 
Calabi-Yau manifold is defined to have no significant warping, such that the warp 
factor is essentially constant. Suppose Tc is the point where the throat is glued to the 
bulk. From Eqs.( ^.l5| ), ( |3.1t)|) and the asymptotic form of the warp factor for large 
r in Eq.( lj.l7j) , we can find an expression for Tc as the size of the throat. Demanding 
h{Tc) ~ 1 gives 


exp(rc/3) 


^/gsM Ue 


(3.20) 


This relation between Tc, e and Is is important to understand the spectrum of the 
graviton in the next sections. 

The physical size of the throat is given by 

I = 6-^A (. 2/3 f ~ \/gsM Is Tc (3.21) 

Jo A(r) 

We assume that I -C , where is the scale of the compactification. This 
can be achieved by a finite value of Tc and Is -C . Note that the warp factor 
notations are related by 


4. KK Spectrum of Graviton in a Warped Deformed Conifold 


As explained in Section 2, the graviton KK spectrum, Eq. (|2.13|) or equivalently 
Eq. (p.l5|) , can be obtained from the following action for a scalar field in ten dimen¬ 
sions 

S = y (4.1) 

where for graviton = Xlm and /r = 0. However, to study the spec¬ 

trum of massive closed strings, we also added the mass term /i which contains the 
contributions of the closed string oscillatory modes to the mass 

= -1)111 (4.2) 

For graviton KK modes, J\f = 1 and g, = 0, as expected. The equation of motion is 

{^/g (m)) - (m) = 0 . (4.3) 

V 9 

We are interested in calculating the KK mass spectrum in the KS throat glued 
to the bulk of Calabi-Yau manifold. To keep the discussion general we consider the 

is transformed to 


metric ( 3.14 ). Using this metric, Eq. 

-^9yi (^y/gd^^ - g^) = 0 


i = 1, ...,6 


(4.4) 
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where is the mass as measured by the four-dimensional observer: (ja) = 


Let us cast Eq.( |4.4|) into a standard Schrodinger equation form that will be 
used to determine the energy spectrum. Upon expanding the Laplace operator in 
Eq.( ^4.4| ) we will get 'L" and where the derivative is with respect to the coordinate 
r. To transform Eg. ( |4.4|) into the form of an ordinary Schrodinger equation, we 
perform a field redefinition like what usually is done in the hydrogen atom. Let us 
decompose the six internal coordinate = 1,...,6 into the “radial” coordinate r 
and hve azimuthal coordinate ^“,0 = 1, ...,5 of the base manifold (like the angular 
coordinates of in Eq. (|3.10|) ) 


^^r.) = A(r)B^^^{T)i>^^y(n ( 4 . 5 ) 

The hnal wave equation in the radial coordinate r will be for function B{t), while 
A{t) is inserted to get rid of the hrst derivative T' term in the wave equation, as 
explained above. The metric in the KS throat takes the following factorized form 


9ab{T, 0) = Kabir) 'fabiO) ■ (4.6) 

where ■jab represents the intrinsic metric component of the azimuthal directions with 
only 9 dependence. Inserting Eq. (^75| , |T^) into Eq. (|4.4|) , we hnd 

^(m)(^) - Veffim, Eyn, t) S(^)(r) = 0 (4.7) 


where 


Veffim, r) = -grr Ejh^/^ " h' + 

\ a,b=l 

and G{t) up to an overall numerical factor is given by 

G{t) ^ . 

The remaining part of the wavefunction is 

a(,t) ^ G(Tr'/\ 


Qab^ 

$ 

^ m 


+ 


Gi/2 


(4,8) 


(4,9) 


( 4 . 10 ) 


The operator is defined by 


^/95 


^ T^ab/—\ 

m,n=l V '5 
5 

a,6=1 


( 4 . 11 ) 
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where ( 75 ) is the (intrinsic) determinant of the azimuthal directions of the throat. 
The operator is closely related to the Laplace operator for the azimuthal direc¬ 
tions. Later on, we shall give more explicitly. Eq.( |4.7| ) is the analog of the usual 
Schrodinger equation with the “energy eigenvalue” £ = 0, while the energy is 
a parameter in the effective potential Veff{n,E^,T). The Schrodinger equation will 
have S = 0 only for a discrete set of values of Em. This results in the quantization 
of the energy E^. Our aim is to solve Eq. ([4.7| ) for a specific geometry of a deformed 
conifold to hnd the KK spectrum, i.e., the discrete values of E^. For some earlier 
works look at | 2 ^. 

We begin with the “S-orbital” solutions of the 10-D graviton, with p = j — 1 = 
Z = 0, In term of the operator this corresponds to 


0 “^$ = 0 . 

(4.12) 

The effective potential is 


,, 1 .4,3 h(T) 

g m + ^1/2 > 

(4.13) 

and G{t), up to an irrelevant normalization, is 


G{t) = (sinh( 2 r) -2r)2/3. 

(4.14) 

The total wavefunction is given by 



(4.15) 


As explained in Section 2, there always exists a zero mode, the 4-D graviton, with 
Tg(r) = constant. This solution corresponds to Bgir) = and Eq. (|4.7| ) is 

trivially satisfied. 


5. Infinite Throat Limit 

For the excited graviton modes the analytical solution of Eq.( |4.7| ) with the above 
effective potential is not known. To get an idea how the spectrum may look like, we 
consider the limit that the throat is extended to inhnity, i.e. Tc —»■ 00 . In this limit 
we can solve the Eq. (|4.7D numerically. 

The spectrum for the hrst hve level are given in the hrst column of table 1. The 
following empirical relation is obtained for the spectrum which is true for all states, 
with an accuracy better than 0.1 percent 

Em = (1.04n+1.01) n = l,2... (5.1) 

where Em = E^gsMl'^e~‘^^^. This indicates that the spectrum are equally spaced, 
a property of the the roots of the Bessel functions. Furthermore, the spectrum is 
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Figure 3: In the left figure Veffir) is plotted iov fj, = j — 1 = I = 0, i.e. excited graviton 
with no angular momentum. In the right figure, two terms of 14// in Eg.( 4.13 ) for the first 
excited state are plotted. It is clear that the last term in 14// quickly reaches the asymp¬ 
totical value 4/9, while the first term quickly falls off. This justifies our approximations 
used in Eq. (|^) . 


quantized in the unit of as one expects in the spirit of RS model, where 

h(0) = Ha is the warp factor. 


We can get useful information by looking at the shape of the effective potential. 
The effective potential is plotted in Figure |^.a. It has a zero hrst derivative and 
a positive second derivative at r = 0. For large r the last term in Eq.( 4.13D dom¬ 
inates and 14// reaches the asymptotical value 4/9. In the Schrodinger equation, 
the effective energy is zero, which means that the “particle” is moving along the r 
axis. The energy E-m appears in 14// and only for some quantized values of the 
wavefunction falls off at large r. 


We can perform some reasonable approximations to obtain an analytical solution 
which captures the above spectrum quite accurately. The two components of 14// 
are plotted in Figure ^.b . The function — J(r)/iF(r)^ begins with a constant value 
~ 0.94 at r = 0 and rapidly fall off as for large r. The function 

begins with the constant value 2/5 and quickly approaches its asymptotic value 4/9. 
It is clear that the second term in 14// dominates for large r and the solution for 
falls off like 
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level n 

conifold 

roots of W 

roots of J 2 

^2 X 53 

1 

2.05 

1.88 

1.71 

2.02 

2 

3.09 

2.98 

2.81 

3.07 

3 

4.12 

4.05 

3.87 

4.11 

4 

5.16 

5.11 

4.43 

5.15 

5 

6.20 

6.16 

5.99 

6.19 


Table 1. In this table for the hrst 5 levels of the graviton 

KK modes with no angular momentum (/i = j — 1 = / = 0) are given. The hrst 
column contains the actual energy spectrum when the metric of the deformed conifold 
Eq. ( |3.14j) is used. In the second column, the approximation Eq. ( |5.2|) is used and E^ 
is determined by roots of Bessel function J^, as in Eq. ( |5.5|) . In the third column, 
the energy is calculated from roots of J 2 . This gives a good order of magnitude 
estimate. In the fourth column, the S'^ x approximation of the deformed conifold, 
Eq. (|5.7|) , is used. It is clear that this approximation is very good. In each case, the 
approximation improves as n increases. 


We introduce the following approximation 
h(r) Xj . 


~ — e 


-At fZv 




4 

9 


6 K{tY 9 ^ GV 2 9 

where Ai = a gsME^\e~‘^l'^ = aE^ and the numerical convention in the hrst term 
are chosen such that the hnal expression takes a simple form. The parameters a and 
V are undetermined yet. They will be determined by the best ht to the spectrum 

Q- 


In this approximation, the diherential equation (^) can be solved analytically. 
The wavefunction, including the pre-factor is 


'h(m) — 


GW 


Xiv 


2T/‘iu 


+ DY, 


Ai z/ 


2t/3.v 


(5.3) 


(sinh( 2 r) — 2 r)h 3 

To have a localized solution for r —>• cx), we obtain W = 0 (since W(^) 
z ^ 0). On the other hand, the denominator factor in 4/(m) vanishes linearly at 
r = 0, corresponding to the fact that the size of 5^ vanishes at r = 0. In order to 
have a regular solution at the bottom of the throat, we obtain, at r = 0 , 

, Ai z/, 




(5.4) 


This indicates that the spectrum is quantized by the roots of the Bessel function. 
The roots Xi^n of W satisfy the following linear relation for large n 


X. 


TITT 


1 TT 

Z/-— . 

^ 2 ' 2 


(5.5) 


18 



























To match with our empirical result Eq.( ^.l|) , this determines 


V ~ 2.44 


a ~ 2.48. 


(5.6) 


This approximate spectrum from Eqs. (|5.4|) , (|5.5|) , ( ^.6|) is given in the second column 
of table 1. The accuracy is very good, and rapidly improves when n increase. The 
reason that v ^ 2 is that the geometry is not pure Ads in the throat. 

This indeed shows that the spectrum is quantized in the unit of / {gsMlD 
like in the RS scenario. In other words, compared to the fundamental scale the 
physical mass is red-shifted by the warp factor /(gsMlg) r\j h{0) = Ha, as 

expected. 

Having obtained the spectrum of the graviton with no angular momentum in the 
azimuthal directions, we like to consider their effects to the graviton KK spectrum. 
One of the technical issues is the geometrical interpretation of the azimuthal direc¬ 
tions of dsg in Eq.( |3.11|) . As mentioned before, at r = 0, g^,g'^ and g^ combine to 
form a round S^, while g^ and g^ form a configuration which is topologically a 5^ 
but shrinks to zero size at r = 0. For large r, on the other hand, these five angular 
coordinates asymptotically matches to a T^’^. To simplify the calculation, we will 
approximate this complicated situation to an easier case when the angular part is 
X S^. We expect this to be a reasonable approximation. The line element of the 
metric has the following approximate form 


ds^ = h ’^rj^udx^dx^ 


^ 4/3 ^ 

-I- h^Kij) 


dr'^ , 2 /'''' 

L37f5(7) + *5' 


dhlo -|- sinh^ 


I' 

2 ' 


dVt\ 


(5.7) 


where now dVt\ and dVL^ represents the line elements on unit 5^ and respectively. 

We can expand the angular dependent part of wavefunction <h(in) in Eq.( ^75l) in 
terms of the spherical harmonics of S‘^ and and respectively. (See 

the Appendix for a review of the S'^ harmonics.) 


^(m) = 


(5.8) 


The operator (9“^, defined in Eq. ([4.11|) , now decompose to products of Laplace op¬ 
erators for and S'^ 




abf 


K{t) I cosh" 


1 I /(/ -h 1) 


sinh"^ 


$ 


(m) 


(5.9) 


It is easy to add the mass term /x in the potential. As explained before, this 
corresponds for higher excited states with oscillation numbers M > 1. The effective 
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potential from Eg . (14.81) becomes 
^4/3 /i(rl 


Veff = 


6 




j'-l 


3K{t)^ cosh2(|) 


/(/ + !) 


3 i^(r)3 sinh (|) 


GV 2 


where now 


G{t) = (sinh(2r) — 2r)’^/® cosh(-) sinh(r) _ 


(5.10) 


(5.11) 


In the limit that cosh(|)ii'(r)3/2 _ ^ _ g above expression for 


effective potential reduces to that of Eq. (|4.13|) . As a check of the accuracy of our 
approximation in taking the deformed conifold to be x S^, we present the first 
few energy levels using the above potential with j — 1 = / = 0 and compare them to 
the previous results in table 1 . 

As before, we need to approximate the above effective potential to obtain analyt¬ 
ical solutions. The functions K{t)~^ cosh(|)“3 and K{t)~^ sinh(|)“3 quickly reach 
the asymptotic value 2. The function /K{tY varies slightly like So it is 

a good approximation to treat these terms as constant in Veff in Eg. (|5.1CI|) . 

To obtain the spectrum, we use the WKB approximation 


d r^/-Veff = (n - ^ + ^)7r 


n = 1,2... 


(5.12) 


where n represents the radial quantization of the energy spectrum in the warped 
throat. When / = 0, Ti = 0 and T 2 is the point where 14 //('r 2 ) = 0. The factor —1/4 
comes from the fact that the wave function vanishes at the center, r = 0. When 
/ 7 ^ 0, the equation Veff = 0 has two roots : ri and T 2 . 

Performing the integral (|5.12|) , we obtain 


E2 — Q 2 _ Q arctan 


2 , =1.03(n-- + ^) (5.13) 


where 


and 


= c d (j^ — 1) -|- / /(/ -|- 1) -|- .65 , 

1 2^3 

/■ = ^5753 ( s , M = — g , MN '. 


(5.14) 


(5.15) 
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n 

h 

J - 1 

1 

Em 

WKB 

1 

0 

0 

0 

2.02 

1.87 

1 

0 

1 

0 

3.61 

3.53 

2 

0 

1 

0 

4.71 

4.76 

1 

0 

2 

0 

5.09 

4.93 

1 

0 

0 

1 

4.00 

3.77 

2 

0 

0 

1 

5.06 

5.10 

1 

1 

0 

0 

3.38 

3.31 

2 

1 

0 

0 

4.57 

4.53 

1 

1.3 

0 

0 

3.97 

3.90 

1 

1 

1 

0 

4.39 

4.30 

1 

1 

1 

1 

5.81 

5.02 


Table 2. In this table, a sample of the states in the energy spectrum is given 
when various quantum numbers are present. Column labeled by contains the 
actual dimensionless value of energy, while the column labeled by WKB contains the 
estimates from the WKB approximation. In general the accuracy of WKB is around 
a few percent. 


A good £t to the data is given by 


c ~ 3.4, d ~ 1.4, /~3.4. 


(5.16) 


In table 2, the result for Em are given. The accuracy of our WKB approximation 
is in general quite good. In the limit that /i = j^ — l = / = 0, Eq.(b.l3|) rapidly 


approaches the linear regime Eq. (|5.1|) . Also, when one of the quantum number p, j or 
I is turned on, Eq. (|5.13|) also reaches a linear regime between E^ and that quantum 
number. 

As an estimate of /i, we use the numerical values applied in KKLMMT, M = 20, 
Qs = 0.1, which give 


fi ~ 1.3 VW - 1. 

So /t is of order unity for the hrst few excited oscillatory modes. 


(5.17) 


6. A Finite Throat Glued to the Bulk 

Now that we have obtained the KK gravity spectrum for an inhnite warped deformed 
conifold, we like to consider the more realistic case when the throat has a hnite size 
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and is glued to the bulk at Tc- We are interested in the magnitude of a typical KK 
mode wave function outside the throat. Towards the end of brane inflation, the KK 
modes are produced in the inflation (i.e., the A) throat. Naturally, these KK modes 
are localized at the bottom of the A throat. Suppose the standard model branes are 
sitting somewhere else, either as D7-branes wrapping 4-cycles in the bulk, or as a 
stack of 3-branes sitting in another throat. In either case, we like to estimate the 
tunneling amplitudes of the KK modes, that is, their amplitudes in the bulk and/or 
in the standard model thraot. If the tunneling amplitude is very much suppressed, 
then these KK modes will mostly decay into gravitons in the A throat, leading to 
an early universe dominated by gravitational radiation, an undesirable situation. In 
this section, we are interested only in order of magnitude estimates. 

The details of the bulk is not important. The compactihed manifold can have a 
number of throats. The key properties of the bulk relevant for our investigation are 
(1) it is either weakly warped or not warped at all; (b) it has a hnite size. Both of 
these features should be quite generic in a realistic model. To simplify the problem, 
we like to consider a toy model which embodies the above features. 

We can obtain a simple compact manifold by gluing an identical conifold to the 
original one at r = T and then perform an orientifold projection. The range of 
coordinate r is 0 < r < T. For 0 < r < Tc, we keep the same warped deformed 
conifold described earlier. Here, at Tc < T, hijc) = 1. For Tc < t < T, we keep the 
condition h(r) = 1, that is, there is no warping. So the region Tc < t <T dehnes our 
bulk. That is, we take the bulk to be part of the original conifold but with h(r) = 1. 
For large t > Tc, r given by Eq. ( p.l 8 | ) is the natural coordinate and the metric is 

~ dx^dx'" + dr'^ -|- ( 6 . 1 ) 


The physical size of the compactihcation is given by 

31/2 


R ~ 


25/6 


(.2/3 gT/S 


( 6 . 2 ) 


Note that, for r > Tc, only the combination is meaningful. For Tc < r < T, 

the effective potential has the same form as in Eq. ([4.13|) , except that h(r) = 1. 
Furthermore, K{t) falls off like 

To further simplify the exercise, we replace v = 2.44 by its nearest integer, 
namely, v = 2 for 0 < r < Tc. This makes the solution more symmetric when the 
throat is glued to the bulk. With this choice of z/, the comparison to our empirical 
result Eq. (|5.1|) sets a = 3. The spectrum for this choice of z/ = 2 and a is given in the 
third column of table 1. The approximation gets better for higher n and is sufficient 
for our purpose here. In addition, we shall consider only S-waves. 

The effective potential now becomes 




e// 




h{T) (GV2)' 


K{t) 


+ 


GI/2 


+ A)/9 0<r<rc 

^ 4^9 Tc<t<T 
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V(T) 



Figure 4: In this figure the qualitative shape of 14 // given in Eq. (|6.3f ) is plotted. 14 // in 
the throat region quickly reaches the vale 4/9. The throat is glued to the bulk at The 
bulk ends at T, the position of the 0-3 plane. There is a mirror image of this potential for 
T <T <2T. 


where A 2 = . The shape of the effective potential is plotted in hg. ^ 

The continnity of the potential at reqnires that 

A 2 = Ai. (6.4) 


This is eqnivalent to Eq.( 3.20 ), reqniring that hir^ = 1. Now, the wave fnnction 
solntion is still the same as in Eq.( p.3|) , bnt with the new bonndary condition at Tc. 
The new wave fnnction becomes 

d'(m) = ('sinp('2r) - 2r)i/3 [<^i,2'>^2(Ai,2e^^/^) + T>i,2E2(Ai,2e^^/3)] (6.5) 


The solntion in the throat with in the arguments of the Bessel functions rep¬ 
resent an exponentially falling off solution, while, in the bulk, the Bessel functions 
with argument represent a radially propagating plane wave, which decays only 
like r“^. This can easily be understood by choosing the more natural coordinate r 
given in Eq. (|3.18|) for the asymptotical behavior of J 2 and l 2 - 

There are 3 length scales in the model, namely, the string scale, the warped 
scale and the size of compactification R. We hnd useful to introduce 3 dimen¬ 
sionless parameters, /3i, {32 and (3^ : 

= , /?2 = Aie"^^/^ , A 2 = ^ = ;323 (6.6) 

Ps 
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where we have also introduced the notation j3ij = (3i/f3j. Here, (3i ~ y/ggMlg ~ 
= Ha is simply the warp factor at the bottom of the throat. (32 ~ 3y/gsM 
which measures the energy (the KK mass) in term of the string scale. Finally 
~ h/R measures the scale of the compactification in terms of the string scale. 
Given /3i and /^s, we like to solve for (32, i.e., the KK spectrum. In order to trust 
the low energy supergravity approximation, -C 1. Also, for small warp factor, 
/3i -C 1. We find that there is a competition between the smallness of (3i and /33. If 
(3i ^ Ps, the spectra are determined by the warped throat; we refer to this as throat 
domination. On the other hand, if /ds -C /di, the spectra are determined by the overall 
size of the compactihcation. This is the large extra dimension limit and we refer to 
it as the bulk domination. In the intermediate region both Pi and Ps play important 
roles. Finally, for the situation of interest, we also have P 2 ^ 1. This is because the 
energy is reduced compared to the natural scale either by warping in the throat 
domination case or by the large volume effect as in the bulk domination case. This 
approximation for P 2 breaks down when the warping or the scale of compactification 
do not create a hierarchy of scales. In this discussion, we do not consider very highly 
excited modes that reach the continuum limit (that is, P 2 ^ 3)- 
Regularity of the solution at r = 0 now requires 


Di = 


J2 (Ai) 

K 2 (Ai) 


Cl 


J2{.P2l) ^ 
Y2{P2i) ' 


(6.7) 


The boundary condition at r = T, corresponding to r = R the position of the 03- 
plane is given by Eq.( p.l4|) . We consider the symmetric solution with respect to the 
03-plane, that is T(m)('r) = T(m)(2T — r). This will set the boundary condition 
T|m)('r = T) = 0, which gives (using the Bessel function recurrence relations) 


D 2 ~ 


■h{.P23) ^ 

Y^{P23) ' 


( 6 . 8 ) 


The continuity of the wavefunction and its derivative at r = Tc will now determine 
the energy, i.e., P 2 , 

„ y3(fe) - n(fe)J2(/32i)i „ , , 

" y2(/32l)|J2(ft)n(/323)-n(/32)J3(/323)| '■ 

and 

J2{P2)Y2{P2i) -Y2{P2)J2{P2i) ^ J2{P2)YPP23) - Y2{P2) UP 23 ) .g ^.3 
Ji{P2)Y2{P2i) — Yi{P2)J2{P2l) jPP2)Y2,{P2p — YPP 2 ) Jli{P23) 

Let us hrst consider the two limits of Eq. (|6.10|) . One limit is the inhnite throat 
studied in the previous section, where the energy spectrum is quantized by the roots 
of J 2 . This corresponds to Tc —>■ T. In this limit /ds —>■ 1 and the denominator of 
right hand side of Eq.(|6.1U|) vanishes. This requires that 


Ji{.P2)Y2{P2i) — Yi{P2)J2{,P2l) — 0 


( 6 . 11 ) 
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Figure 5: In Figure a, (32i (measures the KK mass spectrum) versus /3i3 is plotted by 
solving Eq.( 6.10D numerically. The curve formed by box, big circle and small circle represent 
the first, second and third excited states respectively. Here /3i = 1/50 but the shape of 
the plots do not change if other values for /3i are chosen as input parameter. In the 
asymptotic region when f32i ~ constant, the throat dominates and the spectrum is given 
by Eq.( |6.12| ) (the horizontal lines in b) . The other asymptotic region is when the bulk 
dominates and the spectrum is given by Eq.( 6.14| ) (the curved lines in b). In between the 
spectrum is given by the superposition of these 2 sets of curves. The transition between 
these two regions happens earlier for the excited modes. Also the excited modes halt 
temporarily when f32i approaches the following lower roots of J 2 and afterwards drops as 
given by Eq. (|6.14D . In the right figure the solutions of Eq. (|6.12|) and Eq.( 6.14 ) versus /^is 
are plotted. The agreement between our numerically obtained curve and the superposition 
of the theoretical curves is very good. Note that the energies of different KK modes never 
cross or become degenerate. 


For /92 -C 1 this gives the expected result 


J2(/32i)^0. (6.12) 

The other limit is the no throat case when Tc —> 0, corresponding to a deformed 
conifold without warping. In this limit, /3i —>■ 1 and the left hand side of Eq. (|6.10| ) 
vanishes. This in turn requires that 

J2(/92)l3(/?23) - Y2{I32)UI32^) = 0 , (6.13) 

which for /d 2 <C 1 gives 


W2z) - 0 , 


(6.14) 
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which determines the energy P 2 in terms of the compactification /?3 as expected. 


Condition (|6.13|) is a manifestation of the fact that the wavefunction is smooth at 
r = 0, when the size of vanishes. To summarize, we see that, in terms of /dis, the 
spectrum is given by Eq. (|6.12|) as / 3 i 3 —>■ 0 (the throat domination limit, given by 
the horizontal line in Fig ^b), and by Eq.( |6.14|) as / 3 i 3 —> cx) (the bulk domination 
limit, given by the horizontal line in Fig ^.b). Using Eq. (|5.5| ) for the nth root of J 2 
and and the mth root of J 3 , we find that 


Pi ^ 


4m + 5 
4n + 3 


Ps 5 


(6.15) 


where m and n represents the hrst few excited states. 

The interesting case is the intermediate region between these two limits. Now 
we like to make the following observation: properly interpreted, the spectrum in 
the intermediate region is completely determined by the combination of Eq. (|6.12| ) 
and by Eq.( |6.14|) . To a very good approximation, the spectrum for general Pi^ is 
given by the superposition of the 2 sets of curves shown in Fig ^b. This is checked 
numerically. In Fig.|^.a we have plotted the hrst few KK mass P 21 versus P 13 . In the 
limit that Pi P 3 , the KK mass is quantized by Eq. (|6.12D . As Pi increases and 
becomes comparable to Ps the curve begins to fall down and for Pi 3> P 3 , the KK 
mass is quantized by Eq. (|6.14 ). The agreement between these two graphs is very 
good. 


7. Tunneling and Wavefunction Distributions 

In this section we study the two throat model. The inhation takes place in the A 
throat while the standard model is located in the S throat. A qualitative view of this 
two throat model is sketched in Fig.|^ The effective potential is given in Eq. (|6.3|) . 
Compared to the set up in the previous section we drop the Z 2 symmetry on the 
throats and the throats now have different warp factors. As before, the energy S in 
our Schrodinger equation is zero. There are two barriers given by region ti < t < T2 
and T 3 < T < T/i that the particle must tunnel through. 

The tunneling rate, or transmission coefficient, is straightforward to obtain in 
the WKB approximation. Beginning with the first barrier, the A throat to the bulk, 
the relation between the coefficients of the incoming waves and the outgoing waves 
is given by ||24 



(7.1) 

-«(0A-^) Oa+^J 

= 2exp(^j dr^Veffir)^ 

(7.2) 


26 




















V(T) 



Figure 6: In this figure the qualitative shape of effective potential for the two throats 
model is shown. The inflationary A throat begins at r = 0 and ends at Tc- The standard 
model S throat ends at r^. In between, T2 < t < T3, is the bulk . 


where ri and T 2 are the classical turning point, as shown in Fig. |^. Using Eg. ( |6.3| ), 
we obtain 

e. = 2exp(-2yU^)(|i^3) 

~ . (7.3) 


In the second line the approximation /d 2 1 is taken. The transmission coefficient 

from the A throat to the bulk is given by 


T(A ^ bulk) = 4 ©A + 


0 . 


-2 




(7.4) 


This may be very small. In the throat domination region, using Eqs.( b.lD , (| 
(|OD we obtain 


and 


T( A —>• bulk) ~ {{n + 1) (7.5) 

where Ha is the warp factor in A throat and n is the radial KK mode. 

We are also interested in the probability of tunneling from the A throat to the 
S throat via the bulk, as shown in Fig.|^. The matrix relating the coefficients of the 
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incoming wave from the A throat to the ontgoing wave in the S throat is given by 


1 j 0yi + 0^ i (0A — 0^ ) 


4 —i {Qa ~ ©A^) A ©A^ 


-,-iL 


JL 


©5 + ©-' 
-i (05 - 




©5^ 


Z(05-© 

©5 + ©5^ 


where L is the following phase integral over the bnlk 

r'^3 


L = 


\J-Velf{T)dT 


(7.6) 


(7.7) 


To get some qnantitative idea abont the valne of L, we may assnme that the bnlk is 
symmetric between two throats at r = T, so 

J T2 


\J-Veff{.T) dr 
= 2^/?l3 - 4- 4 arctan(^^^^^SZ^). 


(7.8) 


Also 0A = ©5 (independent of Z 2 symmetry of the bnlk) by noting that 05 is given 
by the same expression as in Eq. (|7.3|) , i.e., 05 is only a fnnction of the KK mass, P 2 , 
which mnst be the same in both throats. This means i7ni(A) = E^i{S), althongh 
the KK qnantization mode ua is in general different from rig. For example when 
both throats are large enongh so Eq. (^T|) holds, we obtain (tta + l)hA = (jis + l)hs, 
where ^5 is the warp factor in the S throat. So the transmission coefficient from the 
A throat to the S throat via a bnlk is given by 

-1 


r( 7 l ^ S) = 4 I (0^05 + COS^ L + (|d + |i )2 L 


-1 


= 4 1^(0^ + cos^ L + 4 sin^ L 
In the absence of the bnlk, L = 0 and T(A —*■ S') is very small. 


T{A ^S)r^ (0a)-" ~ ~ (riAhA) 


A 6 


(7.9) 


(7.10) 


This shonld be compared to the Randall-Snndrum 5-dimensional spacetime scenario 
where T(A —>• S) ~ (nA^-A)" jl^ . 

However, for some KK modes, L = [ni + l/ 2 ) 7 r, so cosL = 0. For these KK 
modes, the transmission coefficient approaches nnity 

T(A^S)~1 (7.11) 

This is the well-known resonance effect. For this resonance effect to take fnll place, 

V® 


2\l P 23 — 4 — 4arctan(- 


-) = [ul + 2 )^ 


(7.12) 
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Suppose we are in the bulk domination scenario and the energy quantization is given 
by Eq. (|6.14|) . For very massive KK mode, /323 ^ 1; using Eq. (|5.5|) for the Bessel 
functions, we hnd that the resonances always occur. For these modes 


n = —1 + 


flL 


(7.13) 


This means that tunneling from the A throat to another throat passing through the 
bulk may not be suppressed at all. In any realistic model, the KK modes would have 
some widths. Suppose hs hA- We expect a u^th KK mode will overlap with a 
small set of u^th KK modes so that the resonance effect will come into play. 

Another quantity of phenomenological interest might be the ratio P of a KK 
mode tE'(m) in the bulk to its value at ri the classical turning point in Fig. |^: 



2 /^2 i 

f J2{(^2)Y2{(32 i) — Y2{(32)J2{P2i)\ 

^(m)(ri) 

16 ' 

[ J2i2)Y2i/32l)-Y2i2)J2i/32l) ) 


(7.14) 


In the throat dominated limit, the term in the bracket scales like and P scales 
like /3| ~ ((n + 1) for the radial excitation n. Since the typical warp factor in 
brane inflation has value ~ 1/5 to 10“^, we see that P can be exponentially small. 
However, in the bulk dominated region, the term in the bracket scales like and 
P increases rapidly and may become comparable to one. We see that there is little 
suppression due to tunneling. 

Another quantity of interest is the ratio of the averaged wavefunction of the 
excited modes compared to the zero mode 




^(m) 


(7.15) 


where 'h{m) is the averaged wavefunction in the bulk. To hnd 77i, we need to use the 
normalized and tl'(m) as given in Eq.(|2.19|). Assuming T Tc 3> 1, one obtains 


f J2{P2)y3{P23) — y2(/?2)»73(/^23) V 
\j2if32)Y2{P2l)-Y2{/32)J2{f32l)J 


In the bulk dominated region, one hnds (using Eq. (|6.13|) ) that 


(7.16) 


77-1 ~ 1 . 


(7.17) 


This is an interesting result, applicable if the lifetime of the KK mode is longer than 
its tunneling rate. It means that, in the large compactihcation limit, this KK mode 
would have a comparable amplitude in the bulk as the zero mode. The reason is 
clear. Although 'h(m) is peaked in the throat, the throat contributes little to the 
volume of the manifold due to the warped measure. As a result of the wave function 
normalization, d'(m) in the bulk should be comparable to that of the graviton. 
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In the throat dominated region, however, one finds 




(7.18) 


which is very small, as one expects. This just indicates that the excited modes are 
confined in the throat. One can also consider the ratio of the wavefunction at r = Tc 
for the excited states compared to the zero mode 


7^2 = 


^(m) ijc) 




(7.19) 


In the throat dominated region 7^2 scales like /Jfg which can be very small for small 
enough Pi. In the bulk dominated region 7^2 reaches the constant value proportional 
to /^Ig, which is the roots of Eq.(p.l4[). In this limit 


7^2 ~ 3, 24,104 


(7.20) 


for n = 1,2 and n = 3 respectively. These numbers are constant and are proportional 


3 ^, where is the n-th root of J 3 from Eq. (|6.14|) . These asymptotic values 


to X 

are independent of choice of Pi, as long as it is in the bulk dominated regime. 


8. Implication for Reheating 

Although how brane inflation takes place in string theory has been studied quite 
extensively recently, the issue of how inflation ends remains a relatively unexplored 
issue (for early analyses, see e.g.. Ref. |^, ^). The collision and annihilation of a pair 
of D3-D3-branes will release energy in the form of massive closed fundamental strings 
as well as D-strings. Besides the very large ones that may end up as cosmic strings, 
most of the energy will end up as low-lying string modes. For closed strings, the 
lowest mode is the graviton. If too much energy ends up in gravitational radiation, 
this will be incompatible with big bang nucleosynthesis. For the brane inflationary 
scenario to be successful, sufficient energy has to be dumped into open string modes, 
i.e., standard model particles. If the standard model branes are sitting at the bottom 
of the inflationary throat, presumably this is not a problem. However, if they are 
sitting in the bulk (say D7-branes wrapping 4 cycles) or at the bottom of another 
throat, then the wave function of the closed string modes must be able to tunnel 
out of the inflationary throat to the standard model throat. This is the heating (or 
pre-heating/reheating) problem in brane inflation. Some recent analyses begin to 
address this issue B i B 0- It is argued that most of the energy will end up in 
massive closed string modes p7|, ^ ^ , which then decay mostly to closed string KK 


modes. It is then important to see how easy it is to transfer the KK energy to open 
strings in another throat. In this paper, we calculate the tunneling rates of S-wave 
KK modes from the inflationary throat to the bulk and to another throat (similar 
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results have been obtained by Chen^). Simple scaling argument shows that they are 
suppressed by higher powers of the warp factor than in the Randall-Sundrum case. 
Very similar properties are expected for the KK modes with angular momenta 
while qualitatively similar properties are expected for the KK modes with S‘^ angular 
momenta. Due to resonance effect, we find that the tunneling rate is very sensitive to 
the details; but, in any realistic scenario, the tunneling of a KK mode from one throat 
to another throat may not be as suppressed as naively thought. This is encouraging 
for a successful exit from brane inflation. However, the result may be very sensitive 
to the specihc properties of the model. Clearly, a more detailed analysis is important, 
the above energy spectrum. 
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A. Spherical Harmonics 


One of the technical issues is the geometrical interpretation of the azimuthal direc¬ 
tions of dsg in Eg ( p.llj) . As mentioned before, at r = 0, g^,g‘^ and g^ combine to 
form a round S^, while g^ and g^ form a configuration which is topologically a 8 “^ 
but shrinks to zero size at r = 0. For large r, on the other hand, these hve an¬ 
gular coordinates asymptotically matches to a T^’^. To simplify the calculation, we 
approximate this situation to an easier case where the angular part is x S^. 

The are harmonics on given by PP 




(A.l) 


where F^\x) ci-re the Fock harmonics, 

F^\x) = (smx)^— 

a(cosx) ^ 

Here, satisfies the following orthogonality condition 


(A.2) 




(A.3) 


^Private discussions. 
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Calling the S'^ harmonics we see that the wavefunction takes the form 

njlml'm' 

Since S'0(4) ~ SO {3) x S'0(3), we may decompose them into Qj „i,m' instead, 
where j = 0,1/2,1,..., and m, m' = —j,..., j. Here j = 2j + 1, with degeneracy 
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